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A Study of the Static-Light B B Parameter 

Joseph Christensen* , Terrence Draper and Craig McNeile^ a * 

a Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506 

We calculate the B B parameter, relevant for B°—B° mixing, from a lattice gauge theory simulation using the 
static approximation for the heavy quark and the Wilson action for the light quark and gauge fields. Improved 
sources, produced by an optimized variational technique, MOST, reduce statistical errors and minimize excited-state 
contamination of the ground-state signal. Renormalization of four-fermion operator coefficients, using the Lepage- 
Mackenzie procedure for estimating typical momentum scales, is linearized to reduce order a 2 s uncertainties. 



1. B B Parameter 

Since the lattice static effective theory has 
fewer symmetries than the full continuum theory, 
when calculating the static-light B B parameter 
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operators besides 0^ must be included. These 
correspond to the following full-theory fermion 
operators (see Flynn et al. H): 
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\ ^fP R q) {h„Pm) 
\ [(jh^P L q) (h^q) + 

+ (bj"P R q) (h„PLq) + 

+2 (bP L q) (bP R q) + 2 (bP R q) (bP L q)] 



- (bP L q) (bP L q) 



(2) 



Os is generated at order a s in the continuum due 
to the mass of the heavy quark. O r and On are 
generated at order a s from the chiral symmetry 
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breaking Wilson mass term. The lattice calcula- 
tion of the static-light B B uses the ratio of two- 
and three-point hadronic correlation functions. 

C 3 (t u t 2 ) N»! 



where the required correlation functions are 
C 3 (h,t 2 ) 

= £ {o\T(x(h,xi)b(0,0)T iq (0,0) 
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The three-point function has a fermion operator 
inserted at the spacetime origin, between two ex- 
ternal B-meson interpolating fields. The times 
are restricted to the range t 2 > > t\. The 
gamma matrices, L} and Tj define the type of 
four fermion operator (see equation ^). A spa- 
tially extended -B-meson operator 



f{r)q(t,x + r)j 5 b(t,x) 



(4) 



is used, where / is a smearing function produced 
by MOST |Q for our static f R study. 

2. Scale Formulation 

Using the integrand of the one-loop perturba- 
tive contribution from the coefficients as a weight- 
ing function, as per Lepage and Mackenzie ||, 
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0.82 



Table 1 

"Typical" operator scales; using (3=6.0 and r=l. 

a "typical" momentum scale can be found (Ta- 
bic |l|). Our value for the scale relevant for (An) 
agrees with that found by Hernandez and Hill H] . 
This is the scale which we claim is relevant for 
this calculation as well. We notice that the scale 
found for Bb is singularly different than the oth- 
ers and claim that each of the other matrix ele- 
ments is describing physics at essentially the same 
scale. However, when a ratio is considered, the in- 
tegrands should cancel, but the scale should not. 
Since the other values are similar, we choose 2.18, 
as it has been used for the f B study. 

3. Calculation of the Coefficients 

The coefficients of the operators are calcu- 
lated by renormalization group techniques. 
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where p , 

a l s , we use a v (q*)=0.18 Q. 

The statistical uncertainties for the coefficients 
are listed in Table ||. There is a systematic error 
due to the linearization of the coefficients which is 
not listed. See reference M for complete details. 



4. Results of Simulation 

The raw lattice B parameters for the operators 
which appear in the lattice-continuum matching 
are determined from a Monte Carlo calculation 
of equation || and listed in Table ||. Table || lists 
the linear combination B n m\ = Bb as a function 

L 

of k and extrapolated to k c using fully-linearized 
tadpole-improved coefficients. For both tables, 
the first errors are statistical (bootstrap) and the 
second are systematic due to choice of fit range. 

We find B B (m b )=0.m±i(3)±l as our calcu- 
lated value. The first two errors are as mentioned 
above. The final error is due to the statistical un- 
certainties in the coefficients. If we run to a scale 
of /i=2GeV, with n/=4, using 

(P0,L- 2 P0,f) 



B B ^) = 



B B (m b ) 



(6) 



^(m b ) 

we find B B (2GeV)=1.05(4). When we convert 
BB(rn b ) to a RG invariant quantity using 

b b = (^H)-( p ^ 2pSj 'b b H (7) 

with 4 flavors, we find B B =1.36(6). With 5 fla- 
vors, we find .Bb=1.40(6). 

5. Comparison to Others 

The simulations using Wilson quarks calculate 
the Bb parameter for quark masses around charm 
and extrapolate up to the physical mass, using a 
fit model of the form B = B° + |£. The value 
of B° ( "extrapolated static" in Table [sj> should 
be the same as the static theory. We scale the 
authors' numbers to 2.0 GcV and 4.33 GeV. The 
JLQCD collaboration cite their A as n/=0 val- 
ues. When Abada et al. quote & Bb for the prop- 
agating Wilson quarks, they use nf=0. When 
we scale these, we list values for both n/=0 and 
n/=4. Soni quotes numbers at 2.0 GeV, but no 
A is given. We use our value for A^ 4 ) as well 
as 200 MeV. We calculate B B (m b ), scale it to 
2.0 GeV using n/=4, and calculate a B B with 
both 4 and 5 flavors. Since all of the "raw" values 
are close to 1.0, differences between the estimates 
of the static Bb are due not to the choice of ac- 
tion, but to choices in the coefficients. See for 
the justification of our choice. 
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Z B L 
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+0.008 
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-0.006 
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+0.008 
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-0.012 
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= -0.137 
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+0.003 






-0.003 


-0.003 


-0.004 
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+0.006 
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-0.005 


-0.007 


-0.008 



Table 2 

The absolute changes from our preferred values of the coefficients Zb l , Zb r , Zb n , and Zb s as the 
parameters q*a, a -1 , m b , ar> d Ac are varied by 10%, and 20%, first individually, and then jointly 
("All"), from our preferred input values. The coefficients are quite insensitive to the particular choice of 
input parameter. 





k = 0.152 


k = 0.154 


k = 0.155 


k = 0.156 


k c = 0.157 


B L 
Br 
Bn 
-IBs 


1.01^(1) 
0.96±i(l) 
0.97±2(3) 
1.00^(2) 


1.02^(1) 
0.96^(2) 
0.961?. (4) 
1.00^(2) 


1.02^(1) 
0.95^(2) 
0.96lJ(4) 
1.0012(3) 


1.031^(2) 
0.95±§(2) 
0.96+^(4) 
1.01±1(3) 


1.031^(2) 
0.95l§(2) 
0.95l 3 2 (5) 
1.01±1(3) 



Table 3 

The raw lattice values for the various Bi parameters: Os has a vacuum-saturation value different from 
that of Ol; with a normalization in which the raw Bi have a common denominator equal to the vacuum- 
saturated value of Ol, (— f-Bs) would identically equal 1.0 if vacuum-saturation were exact. 



k = 0152 k = 154 k = 0.155 k = 0.156 k c = 157 

B B {m b ) 0.95l2(l) 0.961^(2) 0.96+g(2) 0.9811(2) 0.9811(3) 
Table 4 ~ ~ 
Ss(mb) is calculated by combining the raw Bi parameters with the appropriate coefficients. 
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Table 5 



The authors' numbers, quoted at the listed scale, have been scaled to n=2.0GeV and to rrib=4.33GeV. 
If the authors quote a number which we used or reproduced, it is bold-faced in the table. 
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